Cooperative effect of phonons and electronic correlations for superconductivity in 

cobaltates 
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We propose that unconventional superconductivity in hydrated sodium cobaltate Na x CoC>2 re- 
sults from an interplay of electronic correlations and electron-phonon interactions. On the basis 
of the t — V model plus phonons we found evidences for a) unconventional superconductivity, b) 
realistic values of T c and c) the dome shape existing near x ~ 0.35. This picture is obtained for V 
close to the critical Coulomb repulsion V c which separates the uniform Fermi liquid from \/3 x y/3 
CDW ordered phase. 
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The discovery of superconductivity in Na x C 0O2 ■ y H2 O 
at T c ~ 5K for x = 0.35 and y = 1.3l has generated a 
great interest in the solid state physics community. Su- 
perconductivity is observed near doping x ~ 0.35 where 
T c follows a characteristic dome shaped (see also Ref.|jJ). 

Cobaltates may be considered as electron-doped Mott 
insulators with layered structure where the Co atoms are 
in a triangular lattice. Hence, cobaltates are strongly cor- 
related systems and Hubbard ort—J models were derived 
to study these materials^. The importance of strong elec- 
tronic correlation is confirmed by recent photoemission 
studies^ showing for x = 0.3 a reduction of bandwidth 
by a factor of two with respect to the calculated ones& 

Several theoretical works predict unconventional su- 
perconductivity in cobaltates due to strong electronic 
correlationiiL&SiiSiil Recent nuclear quadrupolar res- 
onance (NQR)iS, //SR-i£ and specific healii^ measure- 
ments indicate that superconductivity is unconventional, 
with an order parameter in the spin-triplet channel, ex- 
clude time-reversal symmetry breaking, and are consis- 
tent with nodal lines, such that the recently proposed 
next-nearest neighbor (NNN) /-wave triplet Mate' 1 : " J : 
seems a good candidate. However, estimates of T c lead 
to extremely small values 11 . A natural extension to the 
proposed charge fluctuation mechanism£*il is to consider 
electron-phonon (e-ph) interaction, since charge fluctua- 
tions could be enhanced by it. However, a conventional 
e-ph mechanism would be at odds with unconventional 
superconductivity. 

We show below that the interplay between electronic 
correlations and the e-ph interaction close to a charge- 
density instability lead to i) unconventional pairing, as 
originally proposed«iAi but with ii) a realistic value for 
T c and iii) a dome shape around doping x ~ 0.35. 

We start with the t — V model plus phononsl^: 

H = -t ^ (cUaCja + h.C.) + V Y n i n j 



where t and V are the hopping and Coulomb repulsion 
between the nearest-neighbors (nn) sites i and j on the 
triangular lattice. c\ ia and Ci a are the fermionic creation 
and destruction operators for holes, respectively, under 
the constraint that double occupancy is excluded, rij is 
the fermionic density, af i and aj are the phonon cre- 
ation and destruction operators, respectively. The e-ph 
coupling g as well as the phonon frequency loq are con- 
sidered to be constants. The interaction V is motivated 
by optical absortion experiments^ indicating charge in- 
stabilities at x — 0.25 and x = 0.5, i.e. around the region 
where superconductivity takes place. 

To study the model (1) we will use a recently developed 
path integral large-N approach for Hubbard operatorsil, 
starting with the well known relations^, between the c's 
and the Hubbard operators X; X°° = c] ia and n, = 



Our large-N path integral approach for Hubbard op- 
erators does not require any decoupling scheme as in 
slave bosons, and therefore, problems with gauge fluc- 
tuations or Bose condensation are avoided. At leading 
order, our formalism is equivalent to the slave-boson ap- 
proach. However, at the next to leading order (which is 
necessary to calculate dynamical properties) the two for- 
mulations differ. Our formalism was recently tested by 
comparing dynamical properties with exact diagonaliza- 
tion on small clusters^. 

In order to get a finite theory in the N — > 00 limit, we 
rescale t to t/N, V to V/N and g to g/VN and the spin 
projections a are extended from 2 to N . As in slave- 
boson, to leading order, we have free fermions with an 
electronic band, E k — —tx(cosk x + 2 cos cos \k x ), 
renormalized by Coulomb interactions where x is the 
electron doping. In the following t is considered to be 
1. From this electronic dispersion we obtain a large 
Fermi surface (FS) enclosing the T point. First principles 
calculations 6 predict, apart from this FS, the existence 
of small pockets near the K points. This last picture is 
of interest for some studies were the pockets are relevant 
for superconductivity^. However, it is important to no- 
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FIG. 1: a) Total effective paring V e ff as sum of a pure 
electronic mediated V*7f and a phonon mediated V e 
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teractions. Solid, double dashed, and dotted lines are the 
propagators for fermions, charge fluctuations in the t — V 
model and phonons, respectively. In V*f? , the bare e-ph ver- 
tex, g (solid circle), is renormalized by electronic correlations 
as showed in b). The last diagram contains a four leg vertex 
proportional to g which is generated when our X-operator 
approach is applied to the last term of the Hamiltonian (1). 
The Feynman rules used in the evaluation of the diagrams are 
those presented in Refs. ll7lll9l and they are given in terms of 
X-operators and not in terms of slave particles. 



tice that recent ARPES experiments^ do not show the 
presence of pockets. A theoretical explanation for the 
absence of pockets was given using LDA + U—. 

Using our Feynman rulesiii^ we can calculate the pair- 
ing diagrams of Fig. la. The total pairing effective in- 
teraction Vjff{k — k') is the sum of a pure electronic 



V 



t-V i 



e jj (k — k') and phonon V e ff (k — k 1 ) terms . The 
charge fluctuation propagator (double dashed line) in 
V^ff contains a RPA-type series of electronic bubbles, 
such that VZf is exact in 0(l/N)HLil. 

In Ref.^t|, ^ was shown that the pairing effective in- 
teraction from the pure t — V model is of 0(1/N). Due to 
the rescaling of the e-ph interaction, g, superconductiv- 
ity from phonons also appears at 0(1 /N) and therefore it 
can be treated on an equal footing to superconductivity 
in the pure t — V model. 

Pairing is mediated mainly by charge fluctuations 19 
(double dashed line) in V^jY (first diagram of Fig. la). 
As will be seen later, for the parameters used, the renor- 
malized e-ph coupling (Fig. 3) is still in the weak cou- 
pling regime, so that corrections to the self-energy of the 
phonons can be disregarded as usual. The new contri- 
bution to the e-ph pairing potential is the vertex (dark 
triangle) in the second diagram of Fig. la which repre- 
sents the e-ph interaction renormalized by electronic cor- 
relations. The diagrams of Fig. lb show that the renor- 
malization of the bare vertex g is due to the electronic 
correlations of the pure t — V model, that will give the 
main contribution to our results. 

For a given doping, x, the leading order uniform Fermi 



liquid is unstable against a V3 x -\/3 CDW for V larger 
than a critical value V c . Except for a factor 2, (due to 
different definition of the Coulomb term in the Hamilto- 
nian) our mean-field phase diagram in the V c — x plane 
agrees with the obtained one in Ref. [13 (see Fig. 2 of 
that paper). For instance, in our case, V c ~ 1.1 for the 
doping x = 1/3. However, for x ~ 1/3, the instability 
cannot be seen as a softening of a collective charge mode, 
but a redistribution of spectral weight takes place. For 
the values of V considered here, still a clearly defined 
mode appears at high energy (~ 3t) around the wavevec- 
tor (47r/3,0) but sizeable spectral weight appears on a 
broad structure at low energies with a maximum at a 
scale ~ t/2 (see inset (b) in Fig. 2). The inclusion of 
an antiferromagnetic exchange interaction J brings al- 
most no quantitative changes, as expected for the rather 
hight doping level, and consistent with a paramagnetic 
metal2£. Further details of the density response close to 
the instability will be given elsewhere^. We remind that 
the V c — x line separating the uniform Fermi liquid from 
the v3 x v3 CDW phase has a parabola-like shape with 
a minimum closer to the doping x ~ 0.35 where super- 
conductivity takes place in cobaltates. The proximity 
to a charge instability seems to be confirmed in optical 
experiments^, where an anomalous behavior, such as the 
appearance of a low energy peak is observed. A similar 
situation appears in one quarter filling organic materi- 
als and has been interpreted as being induced by charge 
fluctuations close to the charge-ordered transition^. 
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FIG. 2: (color online) Behavior of V^f f ph {k - k') for differ- 
ent Vs approaching V c . As in Ref. [TUl k' is fixed on the FS 
at the angle m = it/6 and measures the angle of k which 
is always on the FS. When V approaches V c , V^fT becomes 
more anisotropic and more attractive favoring triplet NNNf- 
wave superconductivity. Inset (a): behaviour of V*T/ (k — k'). 
Inset (b): imaginary part of the density response x(Q ) tJ ), a t 
Q = (4/37T, 0), for different V (the same of the main panel) ap- 
proaching V c showing the redistribution of the spectral weight 
discussed in the text. 
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In Fig. 2 we show V*fY and V^ ph for two fc-vectors 

(k and k') on the Fermi surface for x = 1/3 and for 
different Coulomb repulsion V, approaching V c . As in 
Ref. , the momentum fc*' is fixed at the angle 0m = 
7r/6 and k runs around the Fermi surface. The effective 
interaction V*7f v , showed in inset (a), agrees with that 
obtained in Fig. 3b of Ref. |ll|. When V approaches V c , 
the effective interaction V^J^ becomes more anisotropic 
and more attractive, favoring superconductivity in the 
triplet NNN /-wave (see below) in agreement with Ref. 

M 

Figure 2 (main panel) shows our main result. The 
bare e-ph paring effective interaction which is origi- 
nally isotropic (the bare g and luq were considered to 
be constants) becomes strongly /c-dependent due to ver- 
tex corrections by electronic correlations. In the calcu- 
lation of VejTf Ph we used the bare value g obtained from 
A = 2g 2 N(0)/u>Q where N(0) is the bare fermionic den- 
sity of states and, the dimensionless bare e-ph coupling 
A was chosen to be 0.42£. Vertex corrections obtained 
by us using the large-N path integral approach are very 
similar to the early calculation of Ref. [26j which were 
used before for studying transport^i and isotope effect in 
cuprates^. 

The main effect of correlations occurs for V close to 
the critical value V c where the e-ph pairing interactions 
V^jj h is strongly anisotropic following the same shape 
of V^jJ (see Fig. 2). Therefore, charge fluctuations close 
to the V3 x \/3 CDW order enhance the anisotropy of 
V^jT ph with the same symmetry as V^jJ anc ^ reinforce 
it. 

We use the effective potentials to compute the dimen- 
sionless effective couplings in the different pairing chan- 
nels or irreducible representations of the order parame- 
ter. The critical temperature T c can be then estimated 
from T c i = l,13u> c exp(l/Xi), where uj c is a suitable cutoff 
frequency, and i denotes the different pairing channels. 

The effective couplings Aj are defined as 1 ^: 

A = J(dk/\v k \) J(dk'/\v u ,\) 9l (k')V eff (k' -k) gi (k) 

(27T)2/(dk/M)0i(k)a 1 ' 

where the functions <7i(k), encode the different pairing 
symmetries (see Table I of Ref. [ll| f° r the triangular lat- 
tice) , and «k are the quasiparticle velocities at the Fermi 
surface. The integrations are restricted to the Fermi sur- 
face. Aj measures the strength of the interaction between 
electrons at the Fermi surface in a given symmetry chan- 
nel i. If Xi > 0, electrons are repelled. Hence, supercon- 
ductivity is only possible when Xi < 0. 

In Fig. 3 we show Aj as a function of doping for the pure 
t— V-mediated, the renomalized e-ph mediated, and total 
(Xf = A* _y + X e ~ ph ) effective interactions in the most 
relevant channels s-wave and NNNf-w&ve for V = 0.96 
and A = 0.4. For pure electronic pairing, in agreement 
with Ref. ^l]], the lowest effective coupling has triplet 
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FIG. 3: Dimensionless superconducting coupling Aj for the 
pure electronic and phononic models as well as the total case. 
Xi is calculated for the s-wave and NNN f -wave symmetries 
and for V = 0.96 and A = 0.4. The other channels (not 
shown) are not relevant. 

NNN f '-wave symmetry near V c (dashed line in the up- 
per panel). As shown by the dashed line in the upper 
panel of Fig. 3, the lowest value of Aj ~ —0.05 leads to 
exp(— 1/0.05) ~ 10 -9 . This extremely small value means 
that a value for T c of the order of a few Kelvin is only 
possible for an unrealistic high value for the cutoff tu c . 
The situation does not change very much in the most fa- 
vorable case when the system is very close to the y/3x a/3 
CDW order (or V = 1 ~ V c ). Even in this case the value 
of Xi ~ —0.1 is also very unfavorable for superconduc- 
tivity. Experimentally, the dome in T c around x ~ 0.35 
has a finite extent and it is not limited to a single value 
of doping, signaling that the range of interactions V able 
to lead to superconductivity cannot be only limited to 
values very close to V c . 

The situation is different for the phonon-mediated 
case. The solid line in the upper panel of Fig. 3 shows 
A, in the NNNf-w&ve. The absolute values of Xi are 
clearly larger than those from the pure electronic case. 
For x ~ 0.36, Xi ~ —0.22. Furthermore, the phonon- 
mediated case leads to a more pronounced dome shape 
in the doping range around x ~ 0.35. 

In the middle panel we present similar results to those 
of the upper panel for the s-wave symmetry. The solid 
line shows the results for the phonon-mediated case with 
large Aj ~ —0.4. However, the results for the t — V model 
(dashed line) are strongly repulsive, Aj ~ 0.3, diminishing 
the attractive effect given by phonons. 

The lowest panel shows results for Xf . The lowest ef- 
fective coupling is in NNN f -wave symmetry following 
a dome shape around x ~ 0.35 where A^ takes values 
Xf ~ —0.25 which are larger (~ six times) than the cor- 
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responding ones from the pure electronic model. This 
means that T c can be of the order of a few Kelvin for a 
cutoff frequency of the order of a high phonon mode. 
The Debye frequency in cobaltates is expected to be 
large according to recent first principles lattice dynamics 
calculations^ that shows the existence of optical phonons 
as high as 75meV. Certainly it is no our aim to give a 
quantitative value of T c but to show that the present 
treatment leads to realistic scales in critical temperature 
and doping. 

At this point one remark is in order. With increasing A, 
the total superconducting couplings in the s-wave and 
NNN /-wave channels become more attractive and, for 
A > 1 both symmetries are nearly degenerated. However, 
lacking detailed information about the e-ph coupling, we 
take a cautious value for A, that is already sufficient to 
trigger superconductivity, and as shown above, of uncon- 
ventional type. 

As superconductivity has a large contribution from 
phonons we expect a large isotope coefficient. Our theory 
also predicts a rather constant isotope coefficient along 
the dome in contrast to the strong doping dependent iso- 
tope coefficient in cuprates^S. To our knowledge, isotope 
effect experiments are still not available for cobaltates. 

It was recently proposed^!, and studied experimen- 
tally, that hydration causes the electronic structure to 



be more two dimensional. We think that due to this ef- 
fect the Coulomb repulsion V is less screened when the 
system is hydrated. If V is small when the system is not 
hydrated, phonons will favor superconductivity only in 
the s-wave channel. However, the strong repulsion with 
this symmetry from the t — V mediated interaction will 
cancel the pairing from phonons. This may be the reason 
for the nonexistence of superconductivity in nonhydrated 
cobaltate. Note that in Ref. 3] it was shown that T c de- 
creases with decreasing lattice parameter c. Finally we 
would like to stress that the phonon contribution to su- 
perconductivity, differently to the pure t — V case, does 
not need the system to be very close to the charge or- 
dered phase. In fact, we find paring inside the metallic 
phase, in fully agreement with the fact that charge order 
was never observed for doping x ~ 0.35 

In summary, we have shown that a cooperative effect of 
phonons and electronic correlations may lead to uncon- 
ventional superconductivity in cobaltates. We showed 
that superconductivity is possible due to e-ph interac- 
tion where the bare e-ph coupling is renormalized by elec- 
tronic correlations in the proximity of the x -\/3 CDW 
instability of the uniform Fermi liquid. Furthermore, we 
have shown that superconductivity is possible for a real- 
istic value of T c , and we also obtained the characteristic 
dome shape near x ~ 0.35. 
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